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Abstract 

The space quantization induced by a Stern-Gerlach experiment is normally explained by invoking 
the "collapse of the wave function." This is a rather mysterious idea; it would be better to explain 
the Stern-Gerlach results without using it. 

We re-analyze the Stern-Gerlach experiment using path integrals. We find if we model explicitly 
the finite width of the beam, coherent interference within the beam itself provides the space quan- 
tization - without need to invoke the collapse. If we insist on employing only wave functions with 
the space and spin parts kept forcibly disentangled, we recreate the need to invoke the collapse. 

The collapse-free approach makes more specific predictions about the shape and position of the 
scattered beams; if the interaction region has finite length, these may be testable. 

Pending experimental disambiguation, the chief arguments in favor of the collapse-free approach 
are that it is simpler and less mysterious, has no adjustable parameters, and requires the invocation 
of no new forces. 
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INTRODUCTION 



Quantum mechanics is over seventy-five years old but remains something of a mystery. 
Much of the difficulty centers around the location of the Heisenberg cut, the boundary 
between the quantum and classical realms. The main problem is the existence of such a 
boundary in the first place. Given that the world is emphatically quantum mechanical, 
there should be no separate domain of competence for classical physics. All classical results 
should ultimately be explicated in quantum terms. 

The boundary between quantum and classical is normally defined as the point where the 
wave function "collapses." The collapse takes place somewhere between the experiment and 
the observer but is unusual among scientific concepts in that it may be found anywhere but 
where it is being looked for. The doctrine of the collapse of the wave function may be traced 
back to Heisenberg with its mathematical formalization taking place in von Neumann's 
eclated textQ. A rev.ew of the concept and its history is prov.ded by J^flfl. 

While calculations using the collapse have been successful, it creates a number of con- 
ceptual problems |sf| : 

1. The "measurement" problem: We should not have separate physics for measurement 
and process. Measurements are themselves a process and should be described by the 
same physics. Whether we call something a "measurement" or a "process" is a matter 
of labels; it should not affect the physics. 

2. The "event" problem: When does the collapse take place? It must somewhere between 
the quantum and classical parts of the analysis, but where in spacetime is the boundary 
located? If we have N observers each observing the other N — 1 (and themselves for 
good measure) are there N 2 collapses? 

3. The "preferred basis" problem: If we are looking at a collapse along the z axis, then 
how did the system "know" it was to collapse along the z axis? What if it became 
confused and collapsed along the y? 



In 1935 Einstein, Rosen, and Podolsky 6] famously argued that quantum mechanics was 
incomplete in that it did not simultaneously give the values of complementary coordinates, 
say x and p x . Bohr [2]) famously refuted their objection by observing that a theory may 
only be required to predict the results of experimentally measurable quantities (although it 



is free to introduce as many unmeasurable quantities as convenience suggests). Only to the 
extent that x and p x are simultaneously measurable is it required that the theory predict 
them. 

The preferred basis problem represents a different kind of incompleteness. It implies a 
requirement that the quantum mechanic employ intuition in deciding how to analyze a given 
experiment. This implies that quantum mechanics is not formally complete. To put this 
in operational terms: there is not in genera 
the location of the Heisenberg cut. As Pearle 



a way we could program software to define 
^] puts it: the Copenhagen Interpretation "is 



a good practical guide for working physicists. Measurement by an apparatus is like great 
art, we know it when we see it - and we have great artists who can bring it about. But, a 
complete theoretical description? No." 

In the Stern-Gerlach experiment 111 ] the collapse is assigned responsibility for the 



observed space quantization. The Stern-Gerlach experiment is a particularly good place to 
study quantization because in this case many of the usual explanations of quantization are 
not available: 

1. Quantization is often supplied by a container of some kind: a violin in the case of 
musical tones, a "black box" in the case of black body radiation, a microwave cavity 
in the case of some mesoscopic cat experiments, the central atomic potential in the 
case of Bohr-Sommerfeld quantization. But in the case of the classic Stern-Gerlach 
experiment there is no container. 

2. In other cases, quantized inputs naturally produce quantized outputs: if an atomic 
system is going from one quantized energy level to another the radiation it emits in 
so doing is automatically quantized as well. But in the case of the Stern-Gerlach 
experiment the inputs are not quantized. 

3. The program of decoherence jl^l IkI ll^. lisl. llf]] has had significant success in explaining 
the tendency of macroscopic systems to focus on a few pointer states. It has been 
confirmed experimentally, e.g. But decoherence is normally held to occur 



20 



21 



22 1 . In the case of the 



result of interaction with the environment 
Stern-Gerlach experiment we may arrange, in principle at least, for the beam to have 
minimal interaction with the environment. 



3 



In the Stern-Gerlach effect, therefore, we see space quantization without obvious external 
explanation. To reproduce this quantization it would appear we must rely solely on the 
collapse of the wave function, naked and inescapable. 

One reasonable response to the doctrine of the collapse is to assume it is in some sense 
real and to attempt to formalize it well enough to give it "a local habitation and a name." 
This is the approach taken in the spontaneous localization (SL) and continuous spontaneous 
localization (CSL) programs: Ghirardi, Rimini, and Weber 23], Pearled, 24 1, and a good 
general review in Dickson [25]. While these have not been experimentally confirmed 0, 
l26j ]. their analyses have helped to highlight the conceptual difficulties associated with the 
collapse. 



We ourselves were 
tional Interpretation 



ed to consider the problem from the point of view of Cramer's Transac- 
27, 2^. We have found this interpretation fruitful but were troubled by 
the apparent need to invoke the collapse within it, at least to deal with the Stern-Gerlach 
experiment. Given the natural connection between the Transactional Interpretation and 
path integrals we tried using path integrals to analyze the problem. 

We were a bit surprised to find that there was no need to invoke the collapse; space 
quantization appeared spontaneously. Further examination made clear that this result was 
not a side-effect of using the Transactional Interpretation or path integrals; it resulted from 
using wave functions which modeled explicitly the width of the beam and whose spin and 
space components were entangled. 

When, as we will show, we restrict the allowed wave functions to include only those in 
which spin and space components are disentangled, we recreate the necessity to invoke the 
collapse of the wave function. This is reasonable enough. As we know from the program 
of decoherence, the boundary between the quantum and classical realms may be defined by 
the point at which acceptable accuracy may be achieved even if we restrict the solutions to 
those using only disentangled wave functions. 

We will begin by examining [a very simple case of] the Stern-Gerlach experiment in 
classical terms. Then we will employ path integrals to analyze the experiment. We will 
find that space quantization results from coherent interference within the incident beam - 
without need for hand work. Then we will find if we insist the spin and space parts of the 
wave function be disentangled, we reproduce the need for the collapse. With the location 
of the problem tentatively identified as "premature disentanglement," we will then examine 
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FIG. 1: Stern- Gerlach experiment in schematic view 

possibilities for experimental test. These will depend largely on what one means by the 
"collapse of the wave function," a concept which resembles a politician in its ability to 
dodge sharp questions. The most obvious line of experimental attack is simply to note that 
the collapse-free approach makes more specific predictions about the location and shape of 
the resulting wave function. 



CLASSICAL ANALYSIS OF STERN-GERLACH EXPERIMENT 

We begin with the classical problem. This provides us a well-defined starting point, lets 
us establish notation, and gives the value of the action along the classical path. A beam 
of uncharged spin (2m + l)/2 particles is sent in the y direction through an inhomogeneous 

— * — * 

magnetic field B aligned along the z axis. We assume B turns on and off abruptly. Since 

— * 

the Stern-Gerlach effect is produced only by the part of B that varies with z we ignore the 

— * 

constant part of the field. To satisfy the requirement that V-B = the magnetic field should 
include small x and y components, but since we are only interested in points of principle we 
ignore these. 

If the particles have non-zero magnetic moment //, their motion will be governed by the 
Hamiltonian 

H=£-it.H (1) 

and they will experience a force 

F = v(/2-S). (2) 

For definiteness we will take the beam as going in the +y direction. It will appear at 
U = Vai enter the magnetic field at yb, leave at y c , and be detected at yd- We will denote 
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the values of z at y a , y c , and yd as z a (— 0), Zb, z c , and Zd with the corresponding times 
t a = t', tb, t c , and ta- We will assume that the beam is tightly focused and that there are 
no significant p x or p z terms. We will take the magnetic field as going along the z axis and 
having the form 



B (x) pa B z z 
B, i ^ 



dz 



z{0{y-y b ) -0(y-y c ))z 



(3) 



2=0 



where 9 is the unit step function. The non-trivial equations of motion are 



{z,H} 



Pz 



m 



p z = {p z ,H} = fi, 



dB 7 



dz 



(4) 



z=0 



in the interaction region, with fi z the component of the spin along the z axis. The solutions 
are 

dB z 



Pz (t) = Pz (h) + (t~ t b ) Hz 



dz 

2 = 

Z{t)=Zb + {t - tb) — + ^l^z- 
Expressed in terms of the endpoints t b , z b , t c , z c 

t c + 



(t - hf 



(5) 



2=0 



( .s z c - z b dB 
Pz (t) = m- — + n 



t c -t b ' dz 
Az dB, 



= m ^t + ^ dz 



t 

2 = 
(t ~ i) 



2 = 



Zit) = ^ + A^ (t - t6) + 2^^7 



(t - t b ) (t - t c ) 



(6) 



2=0 



To first order in the transit time through the interaction region, the trajectory is a straight 
line through the interaction region followed by a sharp kink up or down, as given by the 
change in p z 



Ap z = P z (t c ) - p z (t b ) = fX; 



dB z 



dz 



At 



2 = 



When p z (tb) is zero the change in z is second order in the transit time 

Az = P -^At + ^At=^^ (At) 2 
m 2m 2m dz n 

2 = 

The velocity is given by 



(7) 



(8) 



v = 



X c f + (Xc - Xbf + (X b - X 1 ) 2 



t-tf 



(9) 




FIG. 2: Classical trajectory in an Stern-Gerlach experiment 



t-t' 
.1 



, Vb-y , 

U = v ^ , t c = V 



At 

y c -y' 



Assuming that the angular deviation induced by the field is small, we have to lowest order 

....y-V a, (iQ) 

(ii) 

(12) 
(13) 



Ay 



Ay . m 
At = — = Ay— 

V P y 

^ A Ay dB z 
=>• Ap z = — fj, z 



z=0 



v dz 

We will usually be taking the limit as Ay — > while the product Ay 

This will result in an angular change of direction of 

Ap z mAy dB z 
A9 = = — — fi z 



2=0 



is constant. 



p y Py dz 
This gives for Zd the position at the detection screen 



(14) 



2=0 



Zd 

L 



LAO 



yd 



y c + yb 



yd-y 



(15) 



In the case of the original Stern-Gerlach experiment the magnetic moment of the silver 
atoms was given almost entirely by the magnetic moment of the unpaired electron 

e _ 



eh _ 

a = -fi b a 

2mc 



fx = s -■ 

mc 
eh 

^ b 2mc 

With a and f3 as the azimuthal and polar angles of the spin vector 



(16) 



z d 



mAy dB 
L — — ji b 



P 2 y 



dz 



2 = 




FIG. 3: Expected versus actual distributions 



mAy dB 
= -L — —li b 



P 2 y 



dz 



cos ((3) 



= -z max cos(/3) 

mAy dB. 
= L — — fi b 



z=0 



Pi 
= v z T 

Ay dB z 

V z = Vb^~ 

Py OZ 

T = L— 

Py 



dz 



z=0 



z=0 



(17) 



where v z is the post-interaction velocity in the z direction and T is the time it takes to 
travel from the interaction region to the detection plane. Isotropically distributed spins are 
distributed as 

— sin (P) dad(3 = - sin (/?) d(3 
Air 2 

= (cos 08)) 

^p(cos(^)) = \ (18) 



giving a flat distribution 



p(z d ) = p(cos(/?)) 



1 



d cos (f3) 
dz d 



2z„ 



(9 (z + 2 max ) -6(Z- ^ m ax)) 



(19) 



And of course what is seen is nothing like that: instead of a flat curve, perhaps gently 
rounded at the ends (if the initial beam is really a Gaussian in z) we see two relatively sharp 
spikes at the far ends of the classical distribution. It is as if the magnetic moment is given 
by [i z = =|=yu b with no intermediate values possible. 



PATH INTEGRAL ANALYSIS OF STERN-GERLACH EXPERIMENT 



We now calculate the results in the Stern-Gerlach experiment without invoking the col- 
lapse. We employ path integrals and limit ourselves to considering only a weak, well-localized 
magnetic field, as above. We will compute the propagator for a general spin-dependent po- 
tential, specialize to the Stern-Gerlach case, then compute the wave function as a function 
of time under the assumption that the initial wave function is a minimum uncertainty Gaus- 
sian. This will be enough to make the qualitative features clear. We will discuss possible 
improvements at the end of this section. 



Path integral calculation of propagator with spin 

We begin with the Schrodinger equation for a particle with spin s represented by a spinor 
with M = (2s + l)/2 components [44], with an M by M Hamiltonian H 

ih^) = H\^) (20) 
=> W{t)) = exp (-^j dt"H{t"))j \^{t')) (21) 
\tp (*)) = \x, s) ( dxdx' Y, & s\ exp (~ ( dt"H (t") \ \x' , s') (f, s' \ if) (?)) (22) 

{«'} V l> / 

where the \x,s) represents some suitable basis in coordinate and spin space. In the case 
we are interested in it may be broken into a kinetic energy part, T, diagonal in the spin 
component and a potential part, V, which in general mixes the spin components 

H = T + V (23) 

For time independent H we may write the propagator K as 

K (t,x,s;t',x',s') = (x,s\exp(-^J dt"H(t") \ \x',s') 

= (x,s\exp{-iH\t-t')/h)\x',s') (24) 

We break the exponential up into a product of exponentials over infinitesimal time intervals 

K (t,x, s;t' , x' , s') = lim (x, s\ (exp (— ieT/%) exp (— ieV /%)) \x' , s') 

N—*oo 



e = tj = t' + je (25) 

The to = f and tw — t times are the starting and ending times respectively. We insert a 
resolution of unity between each term in the product 

K (t, x, s; t', x\ s') (26) 
. N-l 

= lim / dxi . . . dxx-i Yl s j+i\ ex P ( — isT/%) exp (— ieV/h) \xj, Sj) 

{si...sjv-i} 3=0 

Sandwiched between two coordinate states the potential gives 

(x j+1 , s j+1 \ V\xj,Sj) = S (x j+1 - Xj) (s j+1 \ V (xj) \sj) (27) 

We deal with the kinetic energy part of the Green's function by inserting a momentum space 
resolution of unity 

1 = f dpJ2\p,s){p,s\ (28) 
W 

on both sides of the kinetic energy term 

(xj+i, Sj + i\ exp (—ieT/h) \xj, Sj) (29) 

{s} {s>} 



= J d fJ2 j d &J2 (^j+i > s J+i I P> s ) (P> s l ex P f 2^ ) ^' S ' ' Sj/ 



{s} {*'} 

Since 



(p, s | x, s') = , 1 exp ( ^ ) 5 SS ' (30) 
<(2irh) 3 \ h J 



the matrix elements for the kinetic energy are 



(xj + i, Sj + i\ exp (—ieT/h) \xj, Sj) (31) 

Ss i+iaj f fip-x j+1 iep 2 ip-xA 

1 dp exp 1 1 



(2-rch) J \ h 2mh % 

3 / • / - - \ 2N 

m I te ml Xj +1 — Xj \ 



JSj+lSj V 2nieh 6XP [h 2 \ 
We can now write the path integral explicitly 

K (t, x, s; t', x" , s') = 



lim / dxi . . . dxN-i V] 



m 



-3N 



{S1...SJV-1 ) 

xex P (f|(|(^) 2 I-V fi )|| (32) 
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Propagator for Stern-Gerlach experiment 



We specialize this to the case we are interested in: a spin 1/2 particle with a non-zero 
magnetic moment and no charge. The Schrodinger equation is 

d h 2 _ 

ih—ipi(t,x) = -— V 2 ^ (t,x) - ' Bffl^j (*>z) 

{3} 

{j} = {-1/2,1/2} (33) 

Since the magnetic field is pointing along the z axis, we select the representation in which 
the spin is diagonal along the z axis and write this as 

d h 2 - 

iK—if>i(t,£) = -— V 2 ipi (t,x) + J2^B Z (x) {cr z ) ij ^ j (t,x) (34) 

{j} 

We are not insisting that the wave function be either spin up or spin down - that would be 
to invoke the collapse. We continue to see the wave function as including both spin up and 
spin down parts, but treat each part separately for mathematical convenience. The up and 
down components of the wave function will propagate like separate wave functions. At the 
end of the calculation we will join the two halves together. We have 

/exp(- few y gj) ) \ 





ieV\ 


exp y 


h ) 1 



" K J+l JJ \ exp 

(35) 

Since the off-diagonal elements are zero, the sums over the Sj telescope; we may therefore 
write the propagator separately for the spin up and down parts 

K±{t,x;t',x") = (36) 



lim f dx*i . . . dxN-i\hr 
So we may therefore write ip as 



-3N 



exp (f N E (f {^f^f T H b B z (*;))) 



Jdx'K + (t,x;t',x' )ifj + (t',x' \ 
\Jdx'K_(t,x; t',x')iP_ (t',x')J 
Essentially we have taken the coupled equations in ip±i/2 and reduced them to two indepen- 
dent equations. 

Per Schulmanj^. we approximate the propagator as 



K ^~i2^m exp & d<M ' )] (38) 
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FIG. 4: Path integral calculation for Stern-Gerlach 



where S c i(t,t') is the classical action from i! to t and f c i(t,t') satisfies the equations (spe 
cialized from Schulman's treatment to our case) 

d 2 fd(t,t>) n e ., n df d (t,t' 



Tu- 



rn 2 



which implies 



o, f c dt',t') = o, 



f cl (t,t') = t-t' 



dt 



t=t' 



(39) 



(40) 



If the Lagrangian is no worse than quadratic in the integration variables, this is exact. For 
the action we have 



S (t, x; t', x') 



Vl (t) , Pi (*) 

2m 2m 



-dt + 



P 2 At) 

2m 



dt 



+ 



+ 



V 




tc 

I 


(vl (t) 




\ 2m 


tb 




tb 

I 




V 


2m 



9B * m\rn 



(41) 



Substituting in the classical trajectory Eq. (jOJ) we get 

i 1 1 i i — / : ^ IT) 

S (t^ x^ 1 •) x 



m (x — x') m (y — y') ^m {z — z c )" 



2 t - t c 

m (z c - z h f dB z z c + z b ,. 1 
+ —^t. r - ; — + n 



2 t-t' 2 t-t' 

\2 



+ 



2(t c -t b ) r 
m {z b — z'Y 



dz 



(t c - h) 



24m 



/' 



dz 



{t c - hf 



(42) 



2 t b -t' 

We are primarily interested in the value of the z part of the action. Dropping terms 
second order in the magnetic field we have 



m{z — z c ) 2 m(z c - z b ) 



2 t-tr 



+ 



2 t c — t b 



dB 7 , 



dz 



Ay z c + z b ^m(z b - z'Y 



z=0 



v 



2 t b -t> 
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m (z- z c f m (z c - z b f z c + z b m (z b - z'f 

~ 2 t-t c 2 t c -t b TmVz 2 2 t b -t' ^ ' 

Then 



(t, z; t', z') = * I ™ . i / w m / w m x / <fe c d* 6 (44) 

(im(z-z c ) 2 \ (im(z c -z b ) 2 ,mv z z b + z c \ 
x exp — exp — =f i- 



2h t - t c J \ 2h t c - t b h 

( im (z b — z') 2 \ 

x exp 



v 2h t b -t' J 

Recalling t c = t + At/2, t b = i — At/2 and discarding terms of order At, we get 



K^(t,z;t',z') = J .™ A exp I (45) 

/ imv z z(i-f) + z' {t-t) imv 2 z (t - i) (i - t')\ 
x exp Jjf— - _ - — ^— - J 

Again dropping terms second order in the magnetic field 

K (t, x; t>, x') = (t, x; t', x') exp ^rnv z z (t - t')+z^t - t) \ ^ 

In other words the kernel looks as if part of the change in p z takes effect at the beginning 
of the interaction region and the rest at the end (provided the interaction region is not too 
long). If the interaction region is positioned 20% of the way from t' to t, then 80% of the 
change will appear to have taken place at t b , 20% at t c , which makes sense. 

We have taken considerable advantage of the assumed simplicity — weak, constant in 
time, sharply delimited in space — of the magnetic potential. If we were facing a more 
realistic spin-dependent potential, say a time-dependent one, it would make sense to compute 
the path integral using a variational approach. At each time slice we would have a 2 by 2 or M 
by M matrix to diagonalize. The matrix rotations required to do the diagonalizations would 
define a time-dependent "classical" solution between the specified endpoints. Approximating 
the integrals by Gaussians would give the first quantum corrections. Essentially we would 
expect to get the Coherent Internal States solutions of Cruz-Barrios and Gmez-Camacho 



311 ]. if from a slightly different direction. 
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Wave function in Stern-Gerlach experiment 



We take the input wave function as a minimum uncertainty wave packet, a Gaussian with 
initial direction along the y axis at t' 

I / (x ^a) 



(f (f , x' 

ip (if, x' 



exp 



2a 2 



+ ik ■ (x' - x a ) , x a = (0, y a , 0) , k = (0, hk y , 0) 



\<fi(t',x')x-J 
Applying the kernel ()46|) to this and integrating over x! 



(47) 



<p± (t,x) 



{na^fn fit-t') 



exp 



[x x a J 



^ (y-yaY 

2o 2 f(t-t>) ' ■ /(t-tf) , 



+ i 



( 



x exp 



,mi) 



hf(t-f) 



2m f (t - f ) 



fit) 



1 + 



(48) 



This gives a probability of finding a particle of either spin as 
3 



p(t,x) 



exp 



x 2 + (y-(y a + v y (t-t'))Y 



^na 2 \f(t-t')\ 2 v 

x Lp f (*+*.(t-*)r 



exp 



a2|/(t-t')| 2 



|X-I 



l/(*)| S 



1 + 



(49) 



The \f(t)\ 2 terms tells us the smaller the initial width a of the beam the faster it diverges, 
as we would expect from the uncertainty principle. 

In this last equation we see the expected two-humped distribution, with peaks at z — 
=pi> 2 (£ — £). The difference between this and the collapsed version is subtle but definite: 
in the collapsed version the wave function is either spin up or spin down: there is no 
relative phase information between the two. Here the wave function is both spin up and spin 
down. The relative phase information is preserved. There is an exponentially small spin-up 
component to be found in the spin-down direction and vice versa. The collapse-free density 
matrix (looking at the z-dependence only) is given by 

/ \ip + (z)\ 2 \ X+ \ 2 <p* + {z)<p-(z)x* + X-\ 
Ul(*)M*)X-X+ b-W| 2 |X-| 2 / 



^(collapse-free) / \ 



(50) 
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The collapsed wave function has the density matrix 

^-^ w = (i^wi*W ) (51) 

V \ V -(z)\ 2 \ X -\ 2 ) 

The collapse-free and collapsed versions have the same trace and therefore the same proba- 
bility for detection at the detection plane. 

The off-diagonal terms are different. But as <p + is centered on the spin-up beam and 
</>_ is centered on the spin-down beam, the amount of spatial overlap between the two is 
exponentially small. Except near the beginning of the interaction region, the off-diagonal 
terms will be nearly zero. This implies that as a practical matter it would be tricky to 
distinguish between the two. We therefore understand why, even if the collapse-free approach 
is in fact correct, the invocation of the collapse still gives good results; the two density 
matrices are almost identical except near the interaction region. Invoking the collapse of the 
wave function is unnecessary, but almost harmless. The major difference between the two 
approaches is (obviously) that the collapse-free approach does not invoke the collapse. Since 
space quantization still appears, it would seem that the step "collapse the wave function" is 
not necessary to produce that. 

Still, if we have lost a collapse, we have gained a catastrophe. If we monitor p (£, x) = 
\ip (t, x)\ 2 as it moves through the interaction region it will acquire at the beginning a small 
dimple in the middle, as the two spins begin to separate. If we take the centroids of \<p+\ and 
\(p^\ as defining the classical paths, on encountering the magnetic field the initially single 
path splits into two. In the language of catastrophe theory 32| this is a cusp or bifurcation 
catastrophe. 

There is still a sense in which the collapse still takes place. When the particle is ultimately 
detected it will be found definitely at one location, within one of the two Gaussian beams. 
This takes place as far past the interaction region as we have patience to wait for the particle. 
It creates no ambiguity as to what basis to use. 



DISENTANGLING THE COLLAPSE 



If we are correct in thinking that the standard formalism, vigorously applied, will generate 
space quantization without invocation of the collapse, then we would like to understand what 
changes in that formalism might recreate the necessity for the collapse, the better to locate 
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the crux. Our suspicion is that the key problem is the inadvertent projection of classical 
ideas onto the quantum realm. 

It is conventional in discussions of Stern-Gerlach experiments to talk about the spatial 
and spin parts as if they were distinct. The spatial part is treated as a kind of carrier wave; 
the spin as the signal impressed on it. The schematic diagram of the Stern-Gerlach experi- 
ment will almost inevitably look like a road map, with forks at each piece of experimental 
apparatus. It is as if the particles are very small cars, each with a propeller strapped on 
top. Mostly the cars travel along fixed roadways, the pathways determined by the collapse 
of the wave function. At critical junctions, their drivers check which way the propeller is 
pointing at the instant and are so guided, electing the spin up or spin down or the spin 
left or spin right pathway accordingly. While the actual formalism treats spin and spatial 
components as entangled, the more informal discussions that surround each set of formal 
steps often treat spin and space as if they were separate. 

The Schrodinger equation has a (free) spatial part and a spin-mixing part. We will 
therefore attempt to solve this using a deliberately incorrect ansatz (modeled on the lines of 
the intuitive picture). We will write the wave functions as composed of disentangled space 
and spin parts 



In the case of the path integral calculation, the initial wave function had, by assumption, 
this form. But an effect of passing through the magnetic field was to turn it into a more 
complex wave function with different spatial dependence in the spin up and spin down parts. 
With this ansatz, Eq. (J33|) gives 



If the interaction term did not mix the spin and space variables the ansatz would work. If 
the separation of variables were possible, the averages would be the separation constants. 
Using the averages is therefore a reasonable way to handle our incorrect ansatz. If we take 
the average of the spin part as giving the effect of /i on and the average of the magnetic 
field as giving the effect of on \ we get 




(52) 



h—ip(t,x)xi (t) 



2m 



V V (*> x) Xi (t) - fkj -B(x)<p {t, x) Xj (t) 



(53) 




(54) 
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where 

(m (*)> = (xl Ax) = -to Yl Xi (t) OijXj (t) (55) 

(B(t)) = (<p\B\(p) = J dxip j (t, x ) B (x) ip (t, x ) 

The initial value for (ft) is determined by the initial value for the spin, which may point in 
any direction. To lowest order in B, x is a constant, therefore (//(£)) = (//(£')) is, giving 
d h 2 

ih—ip (t, x) = -^VV (*, 5) - (0) • ^ (f ) ^ (t, f) (56) 

At this point the equation for cf> is what we had in the last section for each component of 
if), except that the fx ■ B will span a range of values. Therefore for fixed fl we get the wave 
function in Eq. (1481) with the replacement ~^v z — > — 

'f-fa) 2 .k y (y-y a )\ 



(vra 2 ) 3 / 4 V / (* - *) *V 2a 2 / (t - t>) f (t - t>) ) 

y.™(j { ^ mV * f $-* )z ***** \ (57) 
X6XP vV f(t-t>) 2mf(t-t>) (57) 



with probability 

p(t,x) 



-3 



I / x 2 + {y-{y a + v y {t-t'))) 2 + (z- { -^v z {t-i)) 



\7ta*\f(t-t>)\ 2 a*\f(t)\ 2 

(58) 

This is a Gaussian centered at (t — t), which spans the range from the spin up to 

the spin down values. For isotropic spin distributions it will be distributed as cos{(3). To get 
the full probability distribution we convolute the Gaussian with the flat initial distribution 
for the probability, giving the "expected" classical distribution of Fig. El 

In other words, if we assume the spin and space components are disentangled we repro- 
duce the classical prediction. This makes sense. The program of decoherence suggests the 
boundary between classical and quantum modes of analysis may be defined by the point at 
which it becomes acceptable to replace ent ang led quantum wave functions with disentangled 
classical billiard balls and their cuesticks [45]. 

EXPERIMENTAL PREDICTIONS 

The chief advantage of the collapse-free approach is that it is simpler and less mysteri- 
ous. It appears that if we model the finite width of the beam explicitly and refrain from 
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FIG. 5: Schematic Stern-Gerlach with recombination 



invoking the collapse, but otherwise proceed as usual, we will see space quantization appear 
spontaneously as a result of coherent self-interference effects within the beam. We would 
like to identify experimental differences between the two approaches as well. 

There are two detectible differences between the collapsed and the collapse-free ap- 
proaches. One is that in the collapse-free approach we take the off-diagonal elements of 
the density matrix as non-zero. This implies that the two (or more) separated beams are 
still in principle coherent and if brought together again, would be able to interfere with one 
another. For instance, if we start with a spin +x beam, split it in the z direction, recombine 
the two beams, and check the spin along the x direction, we should find it still in a +x 
state. (See Fig. [S] for a very schematic illustration.) This line of attack as been explored 
by Englert, Scully, and Schwinger in their series of papers 'Is Spin Coherence 

Like Humpty Dumpty?" They show that this is possible, although not particularly easy. 
Maintaining phase coherence while separating and recombining the beams is non-trivial. 

The other (in principle) detectible difference is that the location and shape of the final 
wave function is slightly different in the collapsed and collapse-free cases. The most obvious 
difference comes from the finite length of the interaction region. In the collapsed case the 
beam encounters the magnetic field and at some point collapses and heads either up or 
down, as its spin determines. If the interaction region is Ay long, with a field exerting force 
m ^ z ^§z' on a P a rticle traveling with speed v (see Fig. we will see a change of direction 
given by v z . Per Eq. (jHJ) we have for the displacement in the z direction 

Az = z c - z a = z 



v z . v z Ay 
— At = — — 



(59) 



2 2 v 

And we have that the centroid of the wave function goes in the direction v = (0, v,v z ) 
after this. This lets us backtrack the resultant course, giving an effective location for the 
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FIG. 6: Location of collapse 



point of collapse as 



__^a - A v _ Vb + Vc 

y collapse Vc y c 

v z 2 2 



(60) 



so we have a well-defined location for the collapse, as the center of the interaction region. 
There is no fundamental requirement that y C oiia P se be exactly the average of the start and 
end of the interaction region. If we had been less cavalier in our approximations or looked at 
a more complex situation we would see corresponding adjustments in the value of y C oiiapse- 
The collapse-free approach further predicts that z d , the centroid of the detected beam, will 
be at z d = =f^ (y d - ycollapse)- 

Ironically therefore, it is the collapse-free approach which gives a well-defined location 
for the collapse of the wave function |46j |. It is at the intersection of the backtracked output 
paths, with the paths defined by their centroids. The thus defined collapse is not in general 
actually on either of the paths itself. It is a virtual location (appropriately enough) not 
a real one. In general, in any particular experimental case, the collapse-free approach will 
select one specific location as the point of collapse 4j| • And this location will in turn define 
the expected target points. If the centroids are higher or lower than this, the collapse-free 
hypothesis is falsified. 

The view we are led to by this is that the collapse of the wave function represents 
a kind of approximation scheme. It gives a stick figure representation of a more detailed 
calculation. It cannot itself define the location of the collapse, but it does give a qualitatively 
correct representation. To get a correct analysis the path integrals used above or some other 
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technique must be employed. So long as the width of the beam is modeled explicitly and 
entangled wave functions are allowed, space quantization should be seen. 

We have used the simplest possible treatment consistent with reproducing the qualitative 
features. A number of possible improvements in the treatment of the input wave, magnetic 
field, and detectors are obvious. We could model the input wave not as an elementary 
Gaussian but an arbitrary wave form, perhaps modeled as a sum over Gaussian wavelets 
as described in Kaiser |36j| . We could go to considerably greater lengths in describing the 
magnetic field: use second derivative and higher terms in z (or for that matter, include the 
constant part), include the necessarily present x and y components, turn the field on and off 
a bit less abruptly, allow the field to vary in time, model the magnetic field itself in quantum 
mechanical terms. And we could model the interaction with the detector wave functions 
explicitly, perhaps treating them as [rather narrow] Gaussians. 



DISCUSSION 



We conclude that we do not need to invoke the collapse of the wave function to explain 
the space quantization seen in the Stern-Gerlach experiment. Provided that we model 
the apparatus in sufficient detail, space quantization will appear as a result of coherent 
interference within the beam. 

The term "in sufficient detail" means that we explicitly include the finite width of the 
beam, refrain from limiting the set of allowed wave functions to only those with space and 
spin components disentangled, and model the inhomogeneity in the magnetic field to at 
least first order in z. The minimum width of the beam is of course set by the Heisenberg 
uncertainty principle: if the uncertainty Ap z is to be small enough for a well-defined beam 
to exist, then the width Az has to be at least h/Ap z . The requirement to allow entangled 
wave functions is also obvious. [48] 

We have therefore addressed two of the three problems with the collapse of the wave 
function. Except for considerations of calculational efficiency, there is no preferred basis. 
And we may give a well-defined (if virtual) location to the collapse of the wave function. 
This implies that the formalism - at least for the Stern-Gerlach experiment - is "intuition- 
free" and therefore complete in a formal sense. (We therefore find ourselves in concurrence 

n 

with the position of Fuchs and Peres j23]: "Quantum Theory Needs No 'Interpretation'.") 
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The next question is "If the collapse is not necessary, why has its inclusion not done 
harm?" Our answer is that it represents in most cases a very good approximation. The 
spatial separation of the two spin components will usually take place rapidly. Once this has 
happened, the difference between the calculation with and without collapse is minute. The 
collapse is like homeopathic medicine. It does no good but the difference between using it 
and not using it is so small it does little or no harm. 

Of course homeopathic medicine, while relatively harmless, is not free. The collapse 
complicates analysis both in terms of principle - what exactly is it? - and in terms of practice 
- at what point in the trajectory shall we invoke the collapse? Pending an experimental 
disambiguation, the chief arguments in favor of the collapse-free approach are that is simpler 
and less mysterious, has no adjustable parameters, requires the invocation of no new forces, 
and does not create the "preferred basis" problem. The greater conceptual simplicity of 
the collapse-free approach may be of some benefit in spintronics and quantum computation 
applications. 

We see the collapse of the wave function as defining the boundary between the parts 
of the problem space which must be handled quantum mechanically and those which may 
be treated by classical methods. It is not a physical boundary - indeed the same physical 
region might be treated in some ways quantum mechanically and in others classically (as we 
have done here, by treating the magnetic field in classical terms while treating the particles 
traveling through it quantum mechanically). The collapse is something that exists in our 
approximations, our software, and our minds. It does not, however, exist in the phenomena 
themselves, which are purely quantum mechanical in character. 



We have only dealt with one case of the "co 



looked at the optical Stern-Gerlach effect 
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lapse of the wave function." We have not 
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42| and we have not looked at other 



phenomena involving the collapse of the wave function, i.e. the quantum Zeno effect. It 
would be interesting to see to what extent the approach here may be generalized. We suggest 
as a general rule than whenever there is a force giving a Ap, the corresponding q dimension 
should be modeled explicitly, as a Gaussian distribution or whatever seems appropriate. 

We have also not looked at elaborations on the basic Stern-Gerlach idea. One amusing 
variation is to ask what might happen if the beam is sent through a "sandwich" of alternating 
layers of magnetic fields with slopes pointing in the x and z directions. We note that the 
condition for the beam to split is that the impulse given to the beam ji^-At be significantly 
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greater than the uncertainty in the momentum in the z direction of the beam Ap 



z 



~ jl = h 

' Az a 



42|. If the transition is too smooth, the beam will not be split. But we should suddenly see 



a large number of small beams split off - each new (and sufficiently crisp) boundary should 
produce a new "collapse" - if the condition /i |V-B| ^> -^r is met. An interesting thing here 
is that the condition is dependent on the width of the beam; the wider the beam, the earlier 
we should see the effect. 

While we have focused on the non-relativistic case the same considerations apply to 
quantum electrodynamics. The propagators are different; the principles the same. 
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M. O. Scully, R. Shea, and J. D. McCullen, P 
We are following the treatment in Schulman 
elision of some mathematics. 

In the program of decoherence a second boundary, between reversible and irreversible pro- 
cesses, is crossed at the same point 

It would therefore more accurate to say not that the collapse-free and the collapsed approaches 
give different predictions for y co Uapse but rather that the collapse free approach gives a defi- 
nite prediction and the collapsed approach merely implies that it must be somewhere in the 
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interaction region, in our nomenclature somewhere between yj, and y c . 
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region. This would have precluded attainment of a pure case, off diagonal density matrix, upon 
recombining the beams further down-stream." 43| If achieving an absolute separation between 



the parts of a wave function precluded subsequent interference, the double-slit experiment 
would not work. 
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